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Abstract. We investigate the finite source size effect in the context of the wave optics 
in the gravitational lensing. The magnification of an extended source is presented in an 
analytic manner for the singular isothermal sphere lens model as well as the point mass 
lens model with the use of the thin lens approximation. The condition that the finite 
source size effect becomes substantial is demonstrated. As an application, we discuss 
possible observational consequences of the finite source size effect on astrophysical 
systems. 



1. Introduction 

Gravitational lensing is a characteristic phenomenon of the general relativity and has 
become a very important tool in the fields of cosmology and astrophysics PjJ E|- For 
example, the existence of the massive compact halo objects (MACHO) in the Galaxy was 
revealed by the detection of the lensed amplification of stellar objects and the recent 
measurements of the cosmic shear field provide a constraint on the matter distribution 
independently of the clustering bias [U 03 HE E| ■ Promisingly the gravitational lensing 
will play a more important role with progress in the capability of observational facilities 
in future. For example, it will be a useful probe of the nature of the dark energy 

pa ei maun. 

As a fundamental aspect of the gravitational lensing, the effect of wave optics has 
been investigated by many authors [TJl EH EH E] • Very recently, this subject is revisited 
by several authors, motivated by a possible phenomenon which might be observed in the 
future gravitational wave experiments IIEl HZl HHJ HH1 EOl EU E21 ESI EH ESI EHl EH - 
In the context of the wave optics of gravitational lensing, the argument on the distance- 
redshift relation is also revisited [28] . 

The present paper focuses on the finite source size effect in the wave optics of the 
gravitational lensing. In the first half part of this paper, we present an analytic solution 
for the wave equation with the singular isothermal sphere lens model. In general, it 
is difficult to obtain an exact solution for general lens model, excepting a few simple 
lens models (2HIEH!- Therefore it is useful to obtain such analytic solution for the wave 
equation. The wave optics of the singular isothermal sphere lens has been investigated by 
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Takahashi and Nakamura using a numerical technique [3U1 EH] • We present the analytic 
expression for the amplification factor, which is the first aim of the present paper. Then, 
as an application of the analytic formula, we investigate the finite source size effect in 
the wave optics, which is the other aim of the present paper. Using the analytic formula, 
we consider the energy spectrum from an extended source with a Gaussian distribution 
of surface brightness [22] ■ We investigate the condition that the finite source size effect 
becomes important in the wave optics, including the case near the caustic in the limit 
of the geometrical optics. 

This paper is organized as follows: In section 2, we review the basic formulas for the 
wave optics in gravitational lensing. The limit of the geometrical optics is also reviewed 
for self-containment. Then, we present analytic expressions of the amplification factor 
for the point mass lens model as well as the singular isothermal sphere lens model in 
section 3. In section 4, we investigate the finite source size effect with the use of the 
analytic formulas. In section 5, the validity of the approximation of the point source 
is discussed for the gravitational wave from a compact binary. The femtolensing of the 
gamma ray burst source is also revisited [221 122]; and the finite source size effect is 
considered. The last section is devoted to summary and conclusions. Throughout this 
paper, we use the unit in which the light velocity equals 1. 

2. Review of Basic Formalism 

2.1. Wave Optics under the Thin Lens Approximation 

We consider the background spacetime with the line element, 

ds 2 = g^dafdx" = -(1 + 2U(f))dt 2 + (1 - 2U(f))dr 2 J (1) 

where U(r) is the Newtonian gravitational potential with the condition U(r) <C 1. On 
the Newtonian background spacetime, we consider the wave propagation of the scalar 
field The propagation of the electro-magnetic wave and the gravitational wave can 
be well described by the scalar wave equation [T2*l IT], which is given by 

d^iV^gg^d^) = o. (2) 

This is rewritten as 

(V 2 + c; 2 )$ = Au 2 U(r)<$>, (3) 

on the spacetime with the line element (JJJ), where we assume the monochromatic 
wave with the angular frequency uj. In the present paper we consider the spherically 
symmetric potential. 

It is useful to introduce the amplification factor F (which is called the transmission 
factor in Ref. [JJ) by F = $/$o> where $o is the wave amplitude in the absence 
of the gravitational potential, U — 0. Then, under the thin lens approximation, the 
amplification factor is given by 

A „ (.) poo A 

nu,n) = -r§—£-. d 2 £ex P M(£,^], (4) 
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Figure 1. Configuration of the source, the lens and the observer. cIl, ds and d^s are 
the distances between the lens and the observer, the source and observer, and the lens 
and the source, respectively, fj is the position of the (point) source, and £ is the impact 
parameter. (3 is the unlensed source position angle, 9 is the position angle of the image, 
and the deflection angle is a. This sketch is based on the thin lens approximation that 
the wave is scattered only on the thin lens plane. 



where 4>(£,,fj) is the time delay function (Fermat's potential), which is given by 



d L d s ( i V 
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where di is the distance between the lens and the source, ds is the distance between the 
source and the observer, d^s is the distance between the lens and the source, respectively, 
(see Figure d for the configuration of the lensing system and the definition of variables). 
In general, we may add a term 4> m (v) in the right hand side of © [SO]- However, the 
inclusion does not alter our arguments and we omit it. The two dimensional gravitational 
deflection potential is defined by 

V>(I5=2/ dzU&z). (6) 



Note that \F\ = 1 in the absence of the lens potential U — 0. 

The above formulas can be generalized so as to take the cosmological expansion 
into account. Assuming that the wavelength of the scalar waves is much shorter than 
the horizon scale, Eq. (J3J) is generalized as 

F(uJ,ff) = /j "(I + 'l) r fteKpMl + zS&ff)], (7) 

CbLdLs J-oo 

where di,, ds, and d^s are the angular diameter distances, and Zl is the redshift of the 
lens object. 

It is useful to rewrite the amplification factor F in terms of dimensionless quantities. 
We introduce 

x = f> y = tit^ ^ 



4 



w 



1> 



ds 
dhdhs 
didis 
~dslf 



e (l + z L )u, 



(9) 
(10) 



where £o is the normalization constant of the length in the lens plane, for which we 
adopt 



Co — #E^L, 



where 9e is the Einstein angle, i.e., the solution of the lens equation (fTHjl with (3 — 0. (see 
below.) The effect of the wave optics is characterized by the dimensionless parameter 
w. We also introduce the dimensionless time delay function by 



T(x, y) 



\x-y\ 



ip{x). 



(12) 



Then, the amplification factor is written as 

F ( w ^y) = 7T- d 2 xexp[iwT(x,y)}. (13) 

Z7T2 J-oo 

In the case of the spherically symmetric lens model, the gravitational deflection 
potential if)(x) depends only on x = \x\. Then, the amplification factor is reduced to 
the relatively simple formula 

F(w, y) = —iwe 2 



wy 



dx x J (wxy) exp 



iw 



-x 2 — ip(x) 



(14) 



where Jq(z) is the Bessel function of the zeroth order and y 



2.2. Geometrical Optics Approximation 

In this subsection we consider the limit of the short wave length in the wave optics 
(w ^> 1), which reproduces the conventional geometrical optics in the gravitational 
lensing. In the limit of the geometrical optics, the diffraction integral (|13|) is evaluated 
around the stationary points of the time delay function T(x,y). Thus the stationary 
points are determined by the solution of V x T(x, y) = 0, which is written as 

V = x-V x ip(x). (15) 

This is the lens equation to determine the image position Xj. Eq. (|T3j) is rewritten as 

[3 = 6-3(6), (16) 

where (3 = (Co/^l)^ is the angular position of the source, 6 = {^o/d^x is the angular 
position of the images, and a = (£o/^l) V x ^(x) is the deflection angle (see Figure HJ). 
The time delay function T(x, y) is expressed around the j-th image position Xj as 

T(x,y)=T(x v y) + l £ d a d b T(xj, y)X a X b + 0(X 3 ), (17) 

1 o,6=l,2 

where X = x — Xj. Here, the term in proportion to X vanishes because Xj is the 
stationary point of T(x, y). Inserting Eq. (|T7|) into Eq. (JUJ), we obtain the amplification 
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factor in geometrical optics limit [THl EOl EH] 

F geo (w,y) = lA*(^7')| 1/2ex P iwT(xj,y) 




(18) 
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where the magnification of the j-th image is p(xj) = 1/ det(dy/dxj) and rij = 0,1,2 
when Xj is a minimum, saddle, maximum point of T(x,y), respectively. For the case of 
multi-lensed images, in the geometrical optics approximation, the expression (|18|) means 
that the observed wave is described by a superposition of each wave with the amplitude, 
\fi(xj) | 1//2 , and the phase, wT(xj,y) — (rijir)/2 [T3j . 

3. Analytic Expressions for the Amplification Factor 

In this section we present analytic expressions for the amplification factor (fT^j) for 
the point mass lens model and the singular isothermal sphere (SIS) lens model. The 
expression for the point mass lens model is well known ^21 El El E] ■ Recently the 
amplification factor for the SIS lens model is investigated by Takahashi and Nakamura 
with the use of a numerical method [SDIE]- However, we derive an analytic expression 
for the SIS lens model in the present paper. The SIS lens model is often used to 
study a lensing phenomenon by a galaxy halo and by a cluster of galaxies. In order to 
understand the wave effect in the lensing by a halo, such an analytic formula is useful. 
The amplification factor obtained using the analytic expression for the SIS lens model 
gives us the results consistent with those by Takahashi and Nakamura through their 
numerical method [3TH 13*1]. 

3.1. Gravitational Deflection Potential 

Let us here summarize the relation between the gravitational deflection potential ip{x) 
and the density distribution of a lens object (see e.g., X\). In general the deflection 
potential is given by 



where p(x, z) is the mass density distribution of the lens, which is related to the 
Newtonian potential by V 2 U = AirGp. Thus the lens model is characterized by the 
mass density distribution p(x,z) as well as the gravitational deflection potential ip(x). 

3.2. Point Mass Lens 

The gravitational lensing by a black hole and a compact star is described by the point 
mass lens model, in which we write p{x,z) = M5 <y2 \^)5 <yl \z) 1 where M is the mass of 




(19) 



where S(s) is surface mass density in lens plane, 




(20) 
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Figure 2. Magnification fi(w, y) as a function of dimensionless parameter w for the 
point mass lens model. Here, the source position is fixed as y — 0.1, 0.5, and 1, 
respectively. 



the lens object. Then, the surface mass density is E{x) = M5 (2) (0 = M5 {2) (£ x). In 
this model the characteristic Einstein angle 9e is 



Qe = \ — ;— ; — 3 x 10 D — — — ttt arcsec, (21) 

V d L d s \M Q J V IGpc J y ' 

and the gravitational deflection potential is given by if)(x) = logx. Using a mathematical 
integral formula |3*4] . the expression of the amplification factor (|14|1 yields 

F = e >(^+i° g (™/2)) e > r (l _ ^Ai^l - *-w, 1; -^V), (22) 

where iFi(a, c, z) is the confluent hypergeometric function [HE]- In this model we have 
the dimensionless parameter from Eq. 0, which characterizes the wave optics, 

w = 4GM(1 + z L )u ~ 1.2 x 10~ 4 (1 + z £ ) f — ) . (23) 



,M / UHz, 

Note that w has the meaning of the ratio of the Schwarzschild radius to the wavelength 
of the propagating wave. The wave effect becomes significant when w ~ 0(1). 

We define the magnification by /i(w,y) = \F(w, y)\ 2 , which gives us the expression 



TCW 



1-e" 



i a i 



fi(w,y) = - — l;-wy 2 ) , (24) 



2 



,2 ' ' 2 

from expression (J22|) . The maximum magnification is achieved when y — 0, which 
provides the configuration of the Einstein ring, 

- ™ (o^ 

Umax — ^ _ e _ lru( • l/OJ 

Figure 121 shows the magnification (|2*4^) as a function of the wave characteristic parameter 
w with the source position fixed y = 0.1, 0.5, and 1. For w £ 1, the oscillation feature 
appears due to the interference in the wave effect between the double images (see also 
Figure 0J) . 

We next consider the approximation based on the geometrical optics explained in 
the previous section. The point mass lens model has the two images in the geometrical 
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Figure 3. Magnification as a function of the dimcnsionless parameter w. The solid 
curve is the result of the wave optics fi(w,y), while the dashed curve is /i geo (w,y). 
Here, the source position is fixed as y = 0.5. 



optics. Namely, the lens equation has the two solution (the minimum and the saddle 
points of the time delay function). Then, (|18|) yields 

-1 



F geo (w : y) 





1/2 


exp 




\n- 


1/2 


exp 


iw ^ 



iw(^-(p+ - yf 



log \p+ 
log \p- 
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yVy 2 + 4 
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y Vy 2 + 4 



sin 



v^+4 



y 



VlFTI-y 



(26) 



where the magnification of each image is p,± = 1/2 ± (y 2 + 2)/(2y^/y' 2 + 4) and 
p± = (l/2)(y ± + 4). Then, the corresponding magnification is 

y 2 + 2 



(27) 



Figure El shows the magnification ()24j) and (|2T|) . as a function of the parameter w with 
the source position fixed y = 0.5. For w £ 1, both the curves agree, and the geometrical 
optics is a very good approximation. For w ^ 1, however, the two curves are not in 
good agreement because the geometrical optics approximation is not suitable. 



3.3. Singular Isothermal Sphere Lens 

We next consider the SIS lens model, which can be used for modeling a halo. In this 
model, the density profile is 



p(x, z) 



at 



a; 



2nG(\£\ 2 + z 2 ) 2vrG(|eox| 2 + ^)' 
where a v is the velocity dispersion. Then, the surface density is given by 



a:. 



at 



2G\i\ 2G£ x' 



(28) 



(29) 
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and the gravitational deflection potential is given by ip(x) = x. The Einstein angle of 
the SIS lens model is 

,2 d LS _ w , n - 5 { Ov \ 2 (d LS 



•, = ^ = 5x10-^^1^ (30) 



therefore, the dimensionless parameter w is given by 
w = (1 + z L )cj(4Trcr*) 



2 s 2 dLdLS 



d 



s 



0.01(1 + z L 



o /-i . x/ c» W ^ \fd L dLs/d s \ 

\f_^vYf v \( d L d L s/d s \ , 
j llkm/sJ UhzA lGpc I 



/ ^ \ f d L d L s/d s 
lkm/s7 VlHz/V lGpc 
For the SIS lens model, from Eq. p4jl . the amplification factor is written analytically 
(see Appendix A for derivation), 

F(^) = e'-'i^K'f,F,(l + ; 1; -'-wyA. (32) 



n=0 \ Z Z 

Hence, the magnification is written as 

Mw , B) = |f; m±il ( 2l ^)" /2 iFi /_» 1; < y A |. (33) 

n=o n - v z z / 

The maximum magnification is given by setting y = 0, 

X ' r (l + f) f^^l** 



Y,-^i^(2we^) n/Z \ Z (34) 



n=0 ™ ! 



|l + -(!-«>" W^w 



l + Erf(^(l-<)) 



(35) 



where Erf (2) is the error function (see Appendix A). 

Figure |U shows the magnification fi(w, y) for the SIS model as a function of w with 
the dimensionless source position fixed as y — 0.1, 0.5, and 1, respectively. Note that 
when y — 1 (a single image is formed in the geometrical optics limit), the oscillatory 
behavior appears. Our result is consistent with the previous result [3T?1 l3*Tj. 

Finally in this section let us consider the amplification factor based on the 
geometrical optics estimation. In the SIS lens model, the two stationary points (the 
minimum and the saddle points) appear for y < 1, while only one stationary point 
appears for y > 1. Therefore we have 

f | At+ |l/2 e (-^( 2 /+V2)) _ 1/2^+1/2)) ( y < 1 ) ) 

F geo (w,y) = \ (36) 

{ Kl 1/2 (v > i), 

from Eq. (|18p. where /z± = ±1 + 1/?/. Then, the magnification in the geometrical optics 
is written 



2/y + 2 v /-l + l/y 2 sm(2wy) (y<l), 

Vgeo(u), V) = { (37) 

1 + 1/2/ (2/>l). 
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Figure 4. Same as figure [2] but for the SIS lens model. The behavior is very similar 
to that in Figure |21 




Figure 5. Same as figure |3 but for the SIS lens model. The solid curve is the 
magnification fj,(w,y) and the dashed curve is the corresponding geometrical optics 
formula {j, geo (w,y). Here, the source position is fixed as y = 0.5. The geometrical 
optics approximation is not valid for w ;$ 1. 



Figure 03 compares the magnification and (JHTj) as a function of the parameter w, 
where the source position is fixed as y = 0.5. 

4. Magnification of an Extended Source 

In this section we investigate the finite source size effect in the wave optics. We consider 
the magnification from an extended source with a Gaussian distribution of the surface 
brightness. The analytic formulas for the magnification in the previous section are useful 
in the investigation in this section. 
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Figure 6. Configuration of the gravitational lens system for an extended source. Here, 
d,L, d$ and cIls are the angular diameter distances between the observer and the lens, 
between the observer and the source, between the lens and the sources, respectively. £ 
and if are the dimensional coordinates in the lens and the source planes, respectively. 
rs specifies the position of the source center, and as is the source size. This sketch is 
based on the thin lens approximation. 



4-1. Formulation 



Following Ref . [32 , we consider the integral of the point source magnification weighted 
by the source intensity 

-/ \ S-oo W(y)fj,(w, y)d 2 y 

" Kas ' rs)= n.W(S)<Py ■ (38) 

where we assume the Gaussian distribution of the source intensity 



\ff-Y 



where Y (\Y\ = rs) specifies the dimensionless source position, and as is the 
dimensionless source size. These dimensionless quantities of the source position and 
the source size are related to the dimensional quantities by 

r s = fg a s = &S (40) 
ds^E d s 9 E 

where rs and ds are the source position and the source size, respectively (see Figure 
IHJ). Note that the modified magnification depends on the source size as as well as the 
source position rs- 

We find the magnification can be written as 

r 2 2 

1 T f°° (Ts \ 

H(w,a s ,r s ) = —e 2a s / dyye 2a s J —y fx(w,y), (41) 
a s Jo \a s J 
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where Iq(z) is the modified Bessel function of the zeroth order. 

We will find the analytic expression for the integral (|41|) with the use of the result 
obtained in the previous section. Using the Taylor expansion of the magnification fi(w,y) 
around y = r s , 

CO 1 



n=0 



III 



(42) 



where y) is the n-rank derivative of fi(w,y) with respect to y, the magnification 

([41)1 can be written in the form 



fM(w,a s ,r s )= A n ^ n) (w,y = r s ), 



n=0 



where the coefficient is 

1 



A r . 



a 2 s n\ 



e 2a s 



2 2 

" "™dyye~t%lJ^y)(y-r s ) n . 



(43) 



(44) 



We can evaluate the coefficients in an analytic manner. For example, for the first two 
terms, we have 



A 
A x 



1. 



as 



e 4a s 



1 + 



r s 

a 2 s /r s (r 5 >a 5 ) 
a s ^Jn/2 (r s < a s ) 



24 



4a| 



+ 



2a| 



h 



4a| 



(45) 
(46) 

(47) 



where I\{z) is the modified Bessel function of the first order. The other terms can be 
evaluated in a similar way. Note that the zeroth order term of Eq. ()43|) reproduces the 
magnification of the point source. Hence, our expression (J43)) is based on the expansion 
around the point source limit. 



4-2. Point Mass Lens 

For the point mass lens case, with the use of the expression (J2D), we can evaluate the 
magnification of the extended source. Here, we write the first two terms, 

A i2 {0 \w,r s ) + A l ^ l \w,r s ) 



li(w,a s ,r s ) 



(Jt(w,r s ) 



1 - A x w 2 r s lk 



1^(1 + 2; iwrl 
1; \wr 2 s ) 



(48) 
(49) 



For the point mass lens model, with the use of the approximate expression (|2*7|) . 
the magnification of the extended source is evaluated in the geometrical optics. Figure 
[7| shows p,(w,as,rs) and the corresponding magnification with the geometrical optics 
approximation. Here, in evaluating p,(w,as,rs), we summed the terms up to n = 20, 
and the position of the source center and the radius of the source are fixed as a s = 0.5 
and rs = 0.5, respectively. Note that both the curves agree for w £ 1. Comparing it 
with Figure El the oscillation-amplitude decreases as w becomes large. 
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Figure 7. The magnification of the extended source p(w, as,rs) for the point mass 
lens model, as a function w. Here, the position of the source center and the source 
radius are fixed r$ = 0.5 and as = 0.5. For w £ 1, the wave optics agrees with the 
geometrical optics. 
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Figure 8. The magnification £i(w,as,rs) as a function of w for the point mass 
lens model. Here, the position of the source center is fixed as r$ = 0.5, and the 
three curves show the different source sizes specified by as — 0, as — 0.05 and 
as = 0.1, respectively. As the source size becomes larger, the oscillation-amplitude 
of the magnification decreases. In the computation of this magnification, we used the 
approximation of the geometrical optics. This provides a good approximation as long 
as w £ 1, as demonstrated in figure [7| 

Figure |H1 plots the magnification of the extended source, as a function of w. Here, 
the source position is fixed as r s = 0.5, and the three curves assume the source size 
as = 0, as = 0.05, and a s = 0.1, respectively. As the source size becomes larger, the 
oscillation-amplitude of the magnification decreases. This can be understood as follows. 
The oscillation feature comes from the interference of two waves in the geometrical optics 
in the case of the point source. In the case of an extended source, the wave magnification 
is determined by a superposition of many waves. Then, the clear interference disappears 
by averaging over the phase. 

Now let us evaluate the condition that the finite source size effect becomes 
substantial in an analytic manner. The ratio of the second term to the first term 
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of the right hand side of Eq. (J49|) is 

\F 1 (l + iw, 2; \wr\ 



(50) 



1; {wr 2 s ) 

~ -w 2 r s A 1 + 0(w 4 ). (51) 

The condition that the finite source size effect becomes substantial is \8jl/ n\ ~ 0(1). 
For the case r s 3> 05, we may approximate A\ ~ a 2 s /r s , and we have 

~ -(a s w) 2 . (52) 

Thus, a^w is a key parameter of the finite source size effect in the wave optics. 

Next, let us examine the finite source size effect near the caustic rs = in detail. 
Some aspects have been discussed in Ref. [T5] . The above argument is based on the 
expansion of the magnification in terms of w, which is not suitable for the large value 
of w. We here consider the regime of the geometrical optics. In the limit w ^> 1 and 
y <C l/u> 1//2 , we may write 

fj.(w,y) = nwJ (wy) 2 , (53) 

where we used the mathematical formula [33] 

lim xFi(a, 1; z/a) = h(2y/z) (54) 

a^oo 

with z fixed. Note that the approximate formula ()53)1 is rather general, which can be 
derived from ()14j) with the saddle point method for w ^> 1 ( [13], see also below). 
Substituting (JHH|) into (JUJ, the magnification can be evaluated as 

r| 2 

fL(w,a s ,r s ) = —e 2a s / dy y e 2a s I [—y)j (wy) 2 , (55) 
a s Jo \a s J 

which is valid for rs < as ^ 1/w 1 ^ 2 <C 1. Using the definition of the modified Bessel 
function 

00 z 2m 

(56) 



we can write 



™ 'A „ m (rs/a 2 s ) 2m d m a ((3) 



H(w, a s , r s ) = — e 2a s ^(-1) 



a 2 ^0 H) 2 2 2m dp 

where we defined 



P=l/2a% 



(57) 



ao(0)= [°° dyye-^Mwy) 2 = ±-e- w2 Wl (w 2 /2P). (58) 
Jo 2p 



Using the condition, r s < a s <C l/w 1 ^ 2 <C 1, the first two terms of (jBTjl yield 

I if 
2a 



/*(«;, 05, r 8 ) * 7rwe~ r i/ 2a se- w2a %I (w 2 a 2 s ) f 1 + ^ ) , (59) 
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Figure 9. Magnification p,(w,as,rs) as a function of w. Here, the source size is 
fixed as a s = 0.05, and the position of the source center is rg = 0, 0.025, and 0.05, 
respectively. 
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Figure 10. Magnification p,(w, as, i~s) as a function of w. Here, the position of the 
source center is fixed as r$ — 0, and the source size is as — 0, 0.03, 0.04, and 0.05, 
respectively. 



which reduces to 



H(w,a s ,r s ) 



irwe 



-rf/2a; 




7™ 1 + 0(74/4) 



(was < 1] 



7f 1 

2^7 



-r»/2o| 



1 + 



Iff 

2 4, 



(60) 



^ (i + 0(4/4)) 



a 5 

This expression means the followings: the result of the point source is reproduced for 
was <S 1 and rg = 0, while the finite source size effect becomes substantial for was ^ 1 
and p, approaches to the constant value ^Tr/2/a s for r s = 0. The magnification /2 
becomes smaller as rs/as becomes large. 

Figure M shows fi(w,as,rs) as a function of w with the source size as = 0.05 for 
the different source position rs = 0, 0.025, and 0.05, respectively. In this figure we see 
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that p, approaches a constant value (~ l/a s ) slightly depending on r s , as w becomes 
large. Figure [TD] plots fi(w, as, r s ) with the source position fixed as r s = 0, and the 
different source size as = 0, 0.03, 0.04, and 0.05, respectively. Even for the Einstein 
ring configuration, due to the finite source size effect, the maximum amplification is 
limited by the the factor (~ 1/05) for w ^ l/as- These behaviors are consistent with 
those expected from the above analytic argument. 

The above results demonstrate that the condition that the finite source size effect 
becomes substantial is was > 1. The reason can be understood as follows: As shown 
in Appendix B, the condition, was = 1, is equivalent to the condition that the path 
difference between the PATH A and the PATH B in Figure EH becomes comparable to 
the wavelength. Therefore, the observed wave is a superposition of many waves with 
different phases for was ^ 1- This eliminates the interference feature and decreases the 
oscillation feature in the energy spectrum. The finite source size effect near the caustic 
rs = is also understood in the similar way. The maximum magnification is decreased 
by averaging over the magnification of different phases that depend on the position on 
the source surface. 

4-3. Singular Isothermal Sphere Lens 

We here consider the finite source size effect in the SIS lens model. In this case, the 
magnification of the extended source (JUIJ) can be evaluated with the expression (J3~3~j) . 
Figure ITT1 shows the magnification Ji (solid curve), where we set rs = 0.5 and as = 0.5. 
The dashed curve is the corresponding geometrical optics with ()37|) instead of (}3*3*j) . In 
the numerical computation, we performed the sum with respect to n up to 20. 

Figure H21 plots the magnification of the extended source, as a function of w. Here, 
the source position is fixed as r s = 0.5, and the three curves assume the source size 
as = 0, as = 0.05, and as = 0.1, respectively. As the source size becomes larger, the 
oscillation-amplitude of the magnification decreases. Similarly to the case of the point 
mass lens model, the clear oscillation feature disappears as the source size becomes 
large. In this figure we used the approximation of the geometrical optics in evaluating 
the magnification because of a convenience of numerical technique. The validity of its 
approximation is demonstrated in Figure E] at lease for w ^ 5. The approximation is 
not very good for w ^ 5, however, it does not alter our conclusions. 

We write down the first two terms of the magnification 

jl(w, r s , a s ) ~ A ^°\w, r s ) + A 1 /2 il \w, r s ), (61) 

1 — Aiwrs 

' gglo gggHn) (2we^) n/2 ^(1 - f , 2; jwj) V 

- E~ ^ (2we^/r /2 iiM-f , 1; 1^1) 1 ' 

(62) 
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Figure 11. Same as Figure[71 but for the SIS lens model. For w £ 5, the difference 
between the wave optics and the geometrical optics is negligible. Here, we fixed 
rs = 0.5 and as = 0.5. 
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Figure 12. Same as Figure 03 but for the SIS lens model. Here, the position of the 
source center is fixed rs = 0.5, and the magnification is plotted for 05 = 0, 0.05, and 
0.1, respectively. 



The ratio of the second term to the first term of the right hand side of Eq. (|6*2|) is 

5/2/ \i = — Aitur s 

•g =0 £i ^(-^) (2we^ 2 ) n/2 1 F 1 (1 - I 2; \wr 



x 3? 



1) 



E~ £ ^ 11 (2^/ 2 )"/ 2 1 F 1 (-f, 1; \wr%) 



(63) 



w^VsAj + [2 - - J u^rsAi + 0(^ & / 2 ). (64) 

The condition that the point source approximation breaks is \5ji/ y\ ~ 0(1)- In the 
limit rs ^> a,s, we have Ax ~ Og/rs, then 

" 2,, ,3/2 M fo-l\ n,?n1 _l_ fV„,5/^ 



~ + (2 - 1 w o; + C(^ 5/ )• (65) 

Next, similarly to the point mass lens model, we examine the finite source size effect 
near the caustic rg = 0. Here, let us consider the approximate estimation of Eq. (J 14)) 
using the saddle point method, as demonstrated in Ref. j!5j . Using the approximate 
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Figure 13. Same as Figure El but for the SIS lens model. 
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Figure 14. Same as Figure EH but for the SIS lens model. 



method, for w ^> 1 and !/< 1/ y^, we obtain 

Kw,y)^ . 2lXW n ( * ^ J o(w%*y) 2 , (66) 
|1 - ip"(x*)\ 

where positive solution of the lens equation x = ip'(x). For the SIS model, we 

have 

{i(w , y) ~ 2nw J (wy) 2 . (67) 

Substituting Eq. (|H7j) into (jUJ, we have the same expressions of the magnification as 
Eqs. (J55|) and (|50jl . but with multiplied by the constant factor 2. 

Figure 1121 shows p,(w,as,rs) as a function of w with the source size as = 0.05 for 
the different source position r s = 0, 0.025, and 0.05, respectively, for the SIS lens model. 
Similarly, Figure ITU plots p,(w,as,r s ) with the source position fixed as r s = 0, and the 
different source size as = 0, 0.03, 0.04, and 0.05, respectively. These figures show the 
similar behaviors to those in the point mass lens model. 

Finally in this section, we mention our computation and the numerical convergence. 
We have performed the numerical computation using the package MATHEMATICA. 
The terms in Eq. (J4"3~]) with respect to n are summed up to n = 20. The dashed curves 
in Figure E3 shows each term of A n jj,^ as a function w for the point mass lens model, 
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Figure 15. A Jl /i ( ") in Eq. (|43|l as a function of w for the point mass lens model. The 
dashed curves are the terms for n = 0, 2, 4, 6, 8, 12, 14, and 16, respectively, while 
the solid curve is the sum up to the term for n = 20. Here, we fixed the parameters as 
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Figure 16. Same as Figure El but for the SIS lens model. 



where we adopted the parameters as = 0.05 and rs = 0. The solid curve is the (summed) 
magnification jl(w). Figure IT6l is same as Figure H"5l but for the SIS lens model. As long 
as w £ 20, the convergence of our computation is evident. But for the large value of 
w, our method is not advantageous because higher terms with respect to n is required. 
The numerical methods developed in Refs. [T51 EII1 EI] would be useful for numerical 
computation of general cases. 

5. Discussion 

In this section, let us discuss astrophysical consequences of the result in the previous 
section. We first summarize the conditions that the lensing signature of the wave optics 
may appear in the spectral feature, as follows: 

w ~ 1, (68) 
\5fi/(i\ < 1. (69) 
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The first condition (|68|) is that the wave optics in lensing becomes important for a 
monochromatic wave from a point source. The second condition (|59*jl is that the 
oscillation feature in energy spectra survives against the finite source size effect. Here, 
we consider the case r$ ~ as- Therefore, with Eqs. (p)2*|) and (|65p. combination of both 
the condition simply gives as ~ 1. 

Then, we discuss possible observational consequences of the wave effect in the 
astrophysical situation, the lensing of the gravitational wave from a binary compact 
objects [THj and the femtolensing of the gamma ray burst |321 ESj- 



5.1. Gravitational Wave from a Compact Binary 

For the point mass lens model, the dimensionless parameter w is given by Eq. (J2HJ). The 
condition w ~ 1 yields 

(70) 







( M \ 


(iH Z y 


1 + Zl 


uo 4 m ; 



We also have 

where Hq 1 = (70km/s/Mpc)" 1 = 1.3 x 10 26 m is the Hubble distance. On the other 
hand, for the SIS model, w is given by (J31j) . Then, from w ~ 1, we have 



v 



20 / o v \~ 4 ( d L d LS /d s \ 1 



1Hz/ l + 2 L Vlkm/s/ V Hq 1 



(72) 



We also have 



a s ,6x ID" 11 f (^-] " (?j£) . (73) 

6 V10 3 kmy Vlkm/s/ Vdis/ 

Let us consider the binary compact objects, which consists of the two equal objects with 
the mass m. For the binary we have the relation (Gm) 2 /2L = 4" 2 ^(Gm) 5 / 3 (w/2) 2 / 3 , 
where L is the distance between the two objects and u is the angular frequency of the 
gravitational wave. This equation is rewritten as 

If we set L ~ as, it is clear that as C 1 as long as we consider the gravitational wave 
of the frequency around 1Hz. This means that the point source approximation is very 
good for the gravitational wave of this frequency. 



5.2. Femtolensing 

We next consider the femtolensing, which was pointed out by Gould, and Stanek, 
Paczynski and Goodman [321 dH] • The femtolensing is the lens effect by a tiny mass on 
the gamma ray burst. 
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For the point mass lens model, the condition w ~ 1 is rewritten as 

^^^(JLY 1 (75 ) 
IkeV; l + ^U0 20 g7 ' 1 ; 

where h is the Planck constant. The dimensionless source size may be written as 

/ a s \f M V 1 / 2 / Ho 1 \V2 

This suggests that the finite source size effect is important in the femtolensing by the 
point mass lens. If the source size is larger than 10 5 km, the finite source size effect 
becomes significant. In this case the signature of the interference in energy spectra will 
disappear. 

On the other hand, for the SIS model, (J72|l is rewritten as 

( hu \ 0.8 / a v \-U d L d LS lds V X 

UkeVy l + ^V0.1m/sy V Ho 1 / 

We also have 

b V10 5 W V0.1m/s/ \dr,J K ' 



This suggests that the femtolensing might have occurred due to a very small mass halo, 
if it existed. Such the very small mass halo might be unrealistic in our universe, (cf. 
[TTTj ) . However, Moore et al. have pointed out the possibility of the survival of very 
small mass halos which are produced in the high redshift universe, depending on the 
dark matter model |3H], though the possibility is still open to debate jlOl EH1 Ej • Our 
investigation suggests that the source size is a crucial factor even if the femtolensing 
occurred by such the very small halo. When the source size is larger than 10 5 km, the 
interference signature will be significantly affected by the finite source size effect. 



5.3. Finite Source Size Effect near the Caustic 

We now consider the finite source size effect in the wave optics near the caustic. We 
have demonstrated that is becomes influential when asw £ 1 for as <C 1/y/w <C 1. For 
the gravitational wave from a compact binary, asw <C 1 will be reasonable for general 
situation. Hence, we here consider the femtolensing. We may write 




for the SIS lens model, respectively. These estimations suggest that the finite source 
size effect is important in the femtolensing near the caustic too. 
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6. Summary and Conclusions 

In this paper we investigated the finite source size effect on the wave optics in the 
gravitational lensing. First we presented the analytic expression of the magnification 
for the SIS lens model as well as the point mass lens model. Based on the result, we 
evaluated the magnification of the finite-size source, assuming a Gaussian profile for the 
surface intensity. The analytic expression of the magnification is given in terms of the 
expansion with respect to the source size. This expression is useful to understand how 
the finite source size effect works on the spectral feature of the magnification in the wave 
optics. The condition that the finite source size effect becomes significant is discussed. 
As application of the result, we considered the finite source size effect on the wave optics 
in lensing of the gravitational wave from a compact binary and the femtolensing. For 
the lensing of the gravitational wave, it is demonstrated that the finite source size effect 
can be negligible as long as we consider the gravitational wave of the frequency around 
1Hz. For the femtolensing of the gamma ray burst, we confirmed the result by Stanek 
et al. j^21, for the point mass lens model. We also considered the femtolensing by the 
hypothetically very small halo. The femtolensing might imprint the lensing signature on 
energy spectra if occurred, however, the finite source size effect is crucial. If the source 
size is larger than 10 5 km, the finite source size effect becomes significant and will not 
allow the detection of the interference signature in the energy spectra. 

It is worthy noting the finite source size effect near the caustic. In the wave optics 
of the lens configuration of the Einstein ring, the maximum magnification of the point 
source is not divergent, but is in proportion to the frequency of the wave. But, for 
an extended sources, the finite source size effect becomes substantial for a$w ^ 1, in 
which the maximum magnification is suppressed by the value ^Jir/2/as as long as 
as <C 1/y/w <C 1. This finite source size effect would be influential in the femtolensing 
near the caustic too, if the source size is larger than 10 5 km. 
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Appendix A. Analytic Expression for the Amplification Factor in the SIS 
Model 



In this appendix, we derive an analytic expression of the amplification factor of the SIS 
lens model. Since the gravitational deflection potential is ip(x) = x, Eq. ()14|) is written 

(I o 



F(w,y) 



iwe 



iwy 2 /2 



dx x Jo(wxy) exp 



tw\-x —x 



(A.l) 
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= - iwe lwy2 ' 2 £ /°° cfe J ( W xi/)e^ 2 / 2 . (A.2) 

n=0 n - ^° 

Using a mathematical formula |34j . Eq. ()A.2)) is integrated as 

F(«, y) = £ E<i±£ (2u.e«'^)" /2 ^ (l + " 1; -W) ,(A.3) 

n=0 U - V Z Z / 

where iFi(a,c,z) is the confluent hypergeometric function [S3]. Using the formula 

e-'xFi (-|, 1; *) = (l + |, 1; -z) , (A.4) 

we have 

F( m ,y) = £ Hi±l> (2»e«" 2 )" /2 1; *A. (A.5) 

n=0 n - V Z Z / 

In the case y = 0, the Einstein ring configuration, the amplification factor is 

f°° ( x 2 \ 

?/ = 0) = — iw J dx x exp\iw— — iwxj (A. 6) 

= w J dtt exp f — -t 2 - i 3/2 wtj . (A.7) 



Using a mathematical formula 36J, we have the analytic simple form 

F(w, y = 0) = e~™ /4 D_ 2 (W 4 ^) , (A.8) 

where D^iz) is the parabolic cylinder function. With the use of the error function, 
defined by 

Erf (z) = — f dte~ t2 , (A.9) 

y/TT JO 

we may write 

D- 2 (z) = ^e^z{l - Erf (-|)) - e- 2 / 4 , (A.10) 



and we finally have 



F(w, y = 0) = 1 + -(1 - i)e~ iwl2 ^L 



l + Erf^l-i) 



(a.ii; 



Appendix B. Path Difference and the Finite Source Size Effect 



Here, we consider the path difference between the PATH A and PATH B in Figure 
EJ and show that was = 1 is equivalent to the condition that the path difference is 
comparable to the wavelength. Using the angle of the unlensed source position (3 and 
the angle of the image position 6, the Fermat's potential is given 

mp) = ^(e-p) 2 -m- (b.i) 

In the case f3 = (/3, 0) and 9 = (6, 0), the path difference can be evaluated by 

A0(0, 0) = ***S(0 - 0)(A9 - A(3) - ${0)M, (B.2) 

«L5 
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where we assume A9 <C 9 and Af3 -C ft, and = dip(9)/d9, (see Figure |H1 for the 

definitions of A9 and A/3). 

The gravitational lens equation is given by 



d L d s 



(0-p)-$(0) = 0, (B.3) 
d L s 

and the Einstein angle is defined by the solution of the equation 

d J^Q E _ ^ { e E ) = o. (B.4) 

With the use of the above equations, we have 

Aj>(9,f3) = -ft(9)A(3. (B.5) 
We adopt the approximation for the phase difference, 

A}(9, P)~ - i'(e E )A/3 ~ -^9 E A(3. (B.6) 
Because Af3 = ds/di, the condition |(1 + Zi)A<p\ = X/2tt yields 

w(l + z L )-^as9 E {= a s w) = 1, (B.7) 
d L s 

where A is the wavelength of the propagating wave. 



